Numerical validation of the Kerr metric in Bondi-Sachs form 



Liebrecht R. Venter and Nigel T. Bishop 
Department of Mathematical Sciences, University of South Africa, P. 0. Box 392, Unisa 0003, South Africa 

A metric representing the Kerr geometry has been obtained by Pretorius and Israel. We make 
a coordinate transformation on this metric, thereby bringing it into Bondi-Sachs form. In order to 
validate the metric, we evaluate it numerically on a regular grid of the new coordinates. The Ricci 
tensor is then computed, for different discretizations, and found to be convergent to zero. We also 
investigate the behaviour of the metric near the axis of symmetry and confirm regularity. Finally 
we investigate a Bondi-Sachs representation of the Kerr geometry reported by Fletcher and Lun; 
we confirm numerically that their metric is Ricci flat, but find that it has an irregular behaviour at 
the pole. 
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I. INTRODUCTION 



It is widely believed that the end-state of most massive objects is a black hole described by the Kerr geometry. 
The most powerful source of gravitational radiation is the merger of a black hole with another black hole (or other 
. compact object such as a neutron star). Such processes need to be calculated numerically and the field of numerical 
' relativity has been developed to tackle such problems - see, for example the review Q]. The most popular approach 
r***- . to numerical relativity is based on the ADM formalism, and variations thereof, as discussed in the review cited. In 
■ these formalisms, it is understood how to set initial data for a Kerr black hole 

An alternative approach to numerical relativity is based on an evolution of the metric variables in the Bondi-Sachs 
metric 0, • The approach has been used by a number of different groups - for example 0, IE B 0> IHJ > an d 
for a review see ^lj. Applications involving black holes have been restricted to the Schwarzschild case, because of 
the difficulty in setting initial and boundary data for a rotating black hole. Also, the problem of obtaining just the 
horizon data is much simpler |l3[ . and from such data one could then use a numerical code to construct a complete 
I , Bondi-Sachs representation of the Kerr geometry. 

The representation of the Kerr geometry obtained by Pretorius and Israel (and hereafter this reference will be 
denoted as PI) is not in Bondi-Sachs form. However, the PI coordinates are based on outgoing light cones, and as 
such should provide a useful starting-point for the construction of the Bondi-Sachs form of the Kerr metric. We have 
| constructed a coordinate transformation on the PI coordinates that brings the transformed metric into Bondi-Sachs 
form. In order to validate our metric (and also that of PI) we need to show that the Ricci tensor vanishes. The metric 
is not explicit with respect to the new coordinates, so an analytic approach would seem to be hopeless. Instead we 
wrote code that numerically computes the Ricci tensor for the metric, and confirmed that the metric is Ricci flat. 

In order to describe the Kerr geometry, a metric must not only be Ricci flat but must also be regular on the axis 
of symmetry. Thus, we also find the leading order terms of the metric near the axis, and thereby confirm regularity. 

Fletcher and Lun [T^ have proposed a metric in Bondi-Sachs form that represents the Kerr geometry. We will 
call this the FL metric. The metric has the convenient property that its coefficients are explicit functions of its 
coordinates, so a numerical evaluation of the Ricci tensor is straightforward. This we did, confirming that the Ricci 
tensor vanishes. However, we also found that it does not have a normal behaviour near the pole, thus rendering it 
inappropriate for use in numerical relativity. 

Sec. Illl summarizes results obtained in PI that will be used later. Then in Sec. IIHI we obtain in explicit form the 
equations that relate the Boyer-Lindquist and PI coordinates, and also find explicit expressions for the PI metric 
quantities. Sec. IIVI derives the coordinate transformation needed to change the PI metric into Bondi-Sachs form. 
Sec. discusses the numerical implementation of the computation of the Ricci tensor, and presents the results. The 
properties of the FL metric are discussed in Sec. I VII We end with a Conclusion, Sec. IVIII 

II. BACKGROUND: THE PI METRIC 

The Kerr metric in standard Boyer-Lindquist coordinates (t, r, 6, <fi) is 

,9 S , 9 „ ,„o „9 . 9 „ , ,9 Amar sin 2 9 , , , /_ 2mr \ , 9 

ds 2 = —dr 2 + T,d9 2 + R 2 sin 2 Od(j> 2 d(f>dt - fl — J dt 2 (1) 
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where 

S = r 2 + a 2 cos 2 0, A = r 2 + a 2 - 2mr, R 2 = r 2 + a 2 + 2ma ^ sin2 9 _ (2 ) 

The metric obtained by PI (Eq. (PI-31)) uses coordinates (t, r*, A, (/)), with r* = r*(r, 6*), A = X(r, 9). The metric is 

A L 2 

ds 2 = —{drl - dt 2 ) + —dX 2 + R 2 sin 2 9(d<p - uj B dt) 2 , (3) 
R R z 



where (see Eq. (PI-30)) 



2mar . . 

»B = W , (4) 



and L is defined in Eq. (J5J) below. PI often use instead of A as a coordinate, the two quantities being related by 

A = sin 2 (0») (5) 
(see Eqs. (PI-36) and (PI-37)). Thus an alternative form of the metric © is 

ds 2 = A( rfr 2 _ dt i) + L2 S ™l 29 * d6l + R 2 sin 2 9(tht> - uj B dt) 2 . (6) 
R z R z 

PI define the quantities (see Eqs. (PI-7)) 

P 2 (9, A) = a 2 (A - sin 2 9), Q 2 (r, A) = (r 2 + a 2 ) 2 - a 2 AA, (7) 
and (see Eqs. (PI-36) and (PI-14)) 

f°° dr' f e * dO' 
F(r, 9, A) = / - — / subject to the constraint F = 0, (8) 

Jr Q( r j a) Je "\y > -V 



as well as (see Eqs. (PI-28) and (PI-25)) 



dF 

L = fiPQ with (j, = --^r-. (9) 



Eq. (PI-15), and using Eqs. (PI-9) and (PI-39), defines r* 



• 2 + "\ + r ^ + °'-,? fr '- A v + / jv, A)*- (i») 



A(r) ' r A(r') 

There is one minor difference in the notation used in PI and here. We regard the mass of the black hole m and 
the normalized angular momentum a as parameters: thus we do not show a or m as variables in any functional 
dependence list. 

In order to evaluate the metric, the implicit equations need to be reduced to a numerically solvable form. This is 
now done. 

III. PROPERTIES OF THE COORDINATE TRANSFORMATION AND METRIC VARIABLES 

A. Relation between (r,8) and (r*,0«) 

Eq. © implies 

dr' f 6 * d9' 

0, (11) 



Q(r',X) J e P(9',X) 
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which constitutes a relationship of the form function^, 9, 6*) = 0. We found an exact solution for this equation but 
it was more efficient to solve the r-integral in Eci. (|ll|) numerically and the ^-integral exactly as, 



d0' 



P(9',X) a 
$(a, 0*) — e(cosa)S 

E(z,k) = f 
Jo 



[$(0*,0*)-$(0,0*)] 
„ / sin a 



sm( 
dt 



, sin C/» 



vT^tVi - k 2 t 2 

with the condition that a^O and where e(x) is the signum of x. 
Eq. (|10fl can be written as 

U = h{r)+h{r, A)+/ 3 (M) 
and we now evaluate the integrals I\, li and 1%. 

2m? 

I\ = r + m ln(r 2 + a 2 — 2mr) H i=Qo 



C, 



(12) 
(13) 
(14) 

(15) 
(16) 



where v = m 2 — a 2 ; Qo{x) is the zero-order special Legendre function of the second kind; C is an arbitrary constant 
of integration, and we set C = 0. I2 is expressed as 



h = 



00 >2 



e + a 2 - V(C 2 +a 2 ) 2 ~a 2 AA(C) 



A(C) 



<K 



(17) 



and will be evaluated numerically. I3 reduces to 



h = a 



[E( 7) t)-E(-,t)]-(l-t 2 ) F E ( 7> t)-F E (- )t 



(18) 



. _-, sm9 . . . 

7 = sin — — — , t = sm0* (19) 
sin 

and where £ and _Fg are Legendre elliptic integrals of the first and second kinds (we use the notation Fe, rather than 
the usual F, to avoid any confusion with the F introduced by PI and defined in Eq. JSJl). 

B. The metric variable L 

The function L in the metric needs to be written in a form that can be evaluated explicitly. Now, L = [iPQ and 
the problem is to find an explict expression for /x. Using Eqs. JBJl and © 

^Q(r',A) / , d f sin ' lvI de 

A* = j-i t- i2 = ; 

The first term Ti reduces to 



_ 1 !! v !± ^-w) 

2 7 r [( r ' 2 + a 2 ) 2 - a 2 A(r' 2 + a 2 - 2rnr')] 7 

and this term will need to be evaluated numerically. The second term T2 can be simplified to 

1 a sin0 cos 1 Q(0)e(cos0) - 0(0.) 

T2 = -^^^=^^^= H ;^^^=^^^= (22) 

2 a 2 A(A- l)VA-sin 2 2 a 2 A(A - 1) ^A - sin 2 6» 

where 

n(a) =a[£M)-(l-< 2 ) J F 1 B M)l , r; = sin- 1 ^^, i = sin0*. (23) 
L J sin 0* 
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C. Regularity on the axis 8* = 



The Kerr geometry is axially symmetric, and therefore the metric needs to satisfy a regularity condition on the axis 
6* = 0. In order to check the condition we need to know the values on the axis of the non-zero metric components, 
as well as the leading term in a series expansion in 9* for those components that are zero on the axis. Clearly the 
following metric functions have non-zero limits at 6* = 9 = 0: A,S,i? 2 ,ws; but we will need to work out Q and P 
near the axis so that we can obtain expressions for 9 and L near the axis. 

We assume that 9* and 9 are small so that 



and 



Then Eq. JTTJ leads to 



Next, we evaluate [i as 



and then L = [iPQ is 



P^a^~¥ (24) 



rO, 



y/(r 2 + a 2 ) 



dF _ 1 dF _ r 
"d\ ~ ~WlW* ~ 2a 2 9 2 ' 



(26) 



(27) 



L=?^|±^. (28) 

Thus the non-trivial metric components at, or near, the 0* = axis can be evaluated. We find that the PI metric 
satisfies 

2 2n2 2r a 6lma , nnS 

9e,e„=r, = r 0*, g t <p = - ^ 2 - ^ 2 ■ (29) 

Since 9 2 ge,e, / ' — > 1 as 9* — > 0, the metric is regular on the axis. 

IV. COORDINATE TRANSFORMATION TO BONDI-SACHS FORM 

The transformation 

t -> u = t - r* (30) 
brings the metric |J5J into Bondi-Sachs form except for the last term 

R 2 sin 2 9(dcj)-uJB{du + dr*)) 2 ; (31) 
but this term can be brought into the required form with the transformation 

0-><k = ^ + £r(r«A) (32) 
if the arbitrary function H satisfies the condition 

S -^>=-^. (33) 
or* 
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Substitution of the coordinate transforms (|3U[) and (|32|) into Eq. © results in 

, 2 (-du 2 - 2dudr r )A L 2 sin 2 6 Sf cos 2 6 if d6 2 . , dH(r*, 6U , 2 . 2 

rfs 2 = i — '— + 4 — - + (tu B du + ^ — J -d0* - d(j)*) 2 R 2 sin 2 0. (34) 

ri z it Ou* 

We now need to obtain an expression for dH (£*fi*) _ This j s non-trivial as is given in terms of r and rather than 
in terms of r* and 0* . Since H in the coordinate transform is arbitrary up to the condition i|33|) , the remaining 
component is still undefined. We proceed by defining 



^•• 9 -> = f ATO)" 1 *- (35) 

It then follows that Eq. (|33[) is met, since 



A(s)Q(sA) 
dH(r*,6*) 2mar 



dr* £i? 2 



= -Wb. (36) 



After integration of (J3SJ, \ s obtained as 

where 



WB/3sin(26»«) (37) 



t a to s , . 

/3 = M P 2 + / — da. 38 

Note that the integral in Eq. I|38|) is proportional to T\ defined in Eq. (|21l) and, as remarked earlier, will need to be 
evaluated numerically. 
The metric now becomes 

l2 (-du 2 - 2dudr*)A L 2 sin 2 9*cos 2 9*d9 2 . , .,„, . ,„ /on . 

ds 2 = - -5 h 4 -5 + u B + /3sin(20»)d0» - d<j>*) R sin (39) 

fir fir 

which is the Kerr metric in Bondi-Sachs form. 



A. Regularity on the axis 6* = 

The situation is very similar to that discussed earlier for the PI metric in Sec. IIII CI The only new metric function 
that has been introduced is /3, which evaluates on the axis to 

8r(r 2 +a 2 ) y ' 

The extra terms introduced into the metric components by way of j3 are multipled by 0^ ; an d are therefore too small 
to affect the issue of regularity. 



V. NUMERICAL METHOD AND RESULTS 



A. Coding issues 



We mention only matters that are not routine. The expressions containing elliptic integrals E and Fe were evaluated 
using Bulirsch's elliptical numerical solvers |l5j . The quantity S(z, k) defined in Eq. I|14|) can be expressed as a linear 
combination of elliptic integrals of various kinds, also evaluated using Bulirsch elliptical numerical solvers [l5j . The 
quantity I2 is defined in Eq. (|17|l by means of an integral whose integrand decreases monotonically as £ — > 00, therefore 
standard solvers for infinite boundaries can be used. In similar fashion the first integral in Eq. (jl 1|> and the integral 
in Eq. (|21|l were evaluated using the same Romberg integration methods. 

The quantity /i is singular as 9* — > 0, and this proves to be problematic for accurate numerical evaluation of metric 
quantities near the axis. It is for this reason that, in the Tables below, the smallest value of 9* that is used is 0.3. 



6 



B. Construction of data on a regular grid 



In order to obtain numerically the Ricci tensor of the metric to 0(h 2 ), a 13-point regular grid was used. First we 
find first metric derivatives using the central difference formula, then we calculate the Christoffel symbols, then we 
calculate first derivatives of the Christoffel symbols using the central difference formula, and finally we calculate the 
various components of the Ricci tensor. 

For given (r*,0*) we obtain the Boyer-Lyndquist coordinates (r,8) as follows. We make an estimate for r and use 
Eq. to find 8. Once 8 has been obtained in this manner, (r, 8, 8*) is substituted into Eq. I|15|l establishing a value 
for r*. Of course, in general this will not be the desired value of rv, so we repeat the calculation with a different 
estimate of r, and then use the bisection method to find, to an accuracy of 10~ 12 the value of r which leads to the 
correct value of r*. Now, having obtained (r, 8) for given (r*,0*), we find all the metric coefficients, including L, at 
the point (r*, 



First, the numerical code was validated by using it to find the Ricci tensor of the Kerr metric in standard Boyer- 
Linquist coordinates. Second order convergence to zero was observed. 

The Ricci tensor of the metric <|39[1 was evaluated. Tables H] [H] IIIII and IIVI show values for r, 8, L and (3 on a 
regular (r*,#») grid centered at (0.4,0.3) with m — 1, a — 0.1, and h — Ar» = = 0.01. From this data it is 
straightforward to find numerical values for all components of the metric g a @ on the regular (r* , ) grid, and thus to 
compute the Ricci tensor R a p at the centre of the grid. The results are shown in the second column of Table IVl Also 
in that Table we show R a p evaluated using coarser grids, h — 0.02 and h = 0.04, centered at the same point (0.4, 0.3) 
(but in these cases we do not give the intermediate data r, 8, L or (3). The last columns of the Table show the result 
of convergence testing, and it is clear that all components are second order convergent to zero. 

The values of (r, 8) and of g a p were also found for the cases m = 1, a £ {0.1, 0.2} on regular grids (r», 6**) centered 
at r* £ {0.4, 0.5}5, £ {0.3, 0.6, 0.9, 1.2}, and for three different discretization parameters namely h — 0.01, h = 0.02 
and h = 0.04. Rather than present the convergence rates of all the components of the 16 sets of Ricci tensors, we 
instead calculated the L2-norm defined by 



These results are given in Table EH In all cases, second-order convergence to zero was observed. 

Because the metric (|39f) is derived from the PI metric, the above results imply that the PI metric is also Ricci flat; 
we have performed computations similar to those described above and explicitly confirmed this - see Table IVTTl 



Using coordinates (u, r, 8, 0) Fletcher and Lun [l^] have obtained a metric in Bondi-Sachs form, denoted here 
by FL. The metric coefficients do not contain implicit functions, so a numerical evaluation of the Ricci tensor is 
straightforward, and convergence to zero was confirmed (see Table IVIII(I . However, the metric suffers from an 
important restriction in that it is not regular at 9 = 0. 

The FL metric coefficients at 8 — were found to be 



C. 



Results 




(41) 



VI. THE METRIC OF FLETCHER AND LUN 



/// 



r 2 — a 2 + a 2 tanh 2 a + 2mr 



—r 2 — a 2 + a 2 tanh 2 a 
—r 2 — a 2 + a 2 tanh 2 a 



(r 3 + ra 2 + 2a 2 m) 



ho 




amr tanh 2 a 



-2 
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r 



3 cosh 2 a + 2a 2 m 



(r 2 cosh 2 a + a 2 ) cosh 2 a 



le<f> — 
l<h<h — 



mra 2 sinh 2 a 



(r 2 cosh 2 a + a 2 ) cosh 2 a 

sinh 2 a (r 4 cosh 2 a + r 2 a 2 cosh 2 a + 2mr cosh 2 aa 2 + r 2 a 2 — 2mra 2 + a 4 ) 



(r 2 cosh 2 a + a 2 ) cosh 2 a 

(42) 



where a = a(r, a, to) is dchned by 



dv 



yV 4 + a 2 v 2 + 2a 2 mv 



(43) 



In general, all the angular components of the FL metric are non-zero at 6 = 0, so the metric cannot be regular 
there. Now, ^1 suggests that the axis of symmetry is not at 9 = but its location is a function of r, explicitly 
9 = arcsin(— tanh(a)). We do not comment on the analytic validity of such a statement. Our aim is to find a form 
for the Kerr metric that can be used in numerical codes based on the Bondi-Sachs metric, and whatever the case - 
irregularity at 9 = 0, or an axis of symmetry that is not a straight coordinate line - the metric of [L2T ] is unsuitable 
for this purpose. 



VII. CONCLUSION 



We have computed numerically the Ricci tensors of the metrics (|39|) . PI and FL. In all cases, second order conver- 
gence to zero was observed. We also investigated regularity near the axis = (or 9 — in the FL case), finding 
that the metrics (13911 and PI were regular, but that FL was not. 

The aim of this work is to find a metric in Bondi-Sachs form that represents the Kerr geometry and that can be 
used in a numerical computation. The metric (|39[) partially meets that goal. In order to meet it completely, the 
radial coordinate r* will need to be transformed to an area coordinate say r** (i.e., so that the area of the 2-surface 
u = r** = constant is 47rr 2 l J; the procedure for doing so has been described in Also, there will need to be a 

transformation to the angular coordinates (stereographic) actually used in the codes. Finally, the numerical difficulty 
of finding the metric coefficients near the axis will need to be overcome. While not trivial, this is all in principle 
straightforward, and is deferred to further work. 
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TABLE I: r on a regular (r, = 0.4 + ih, 9* = 0.3 + jh) grid with h = 0.01, m = 1 and a = 0.1; used in the metrics and PI. 



i -2 -1 1 2 

3 

-2 2.3701116298499 

-1 2.3685360418458 2.3701098379124 2.3716892173938 

2.3669659799190 2.3685341996528 2.3701079904032 2.3716873645652 2.3732723345139 

1 2.3685323027830 2.3701060880593 2.3716854567442 

2 2.3701041316400 



TABLE II: 6 on a regular (r» = 0A + ih,6 t = 0.3 + jh) grid with h = 0.01, m= 1 and a = 0.1; used in the metrics and PI. 



i -2 -1 1 2 

3 

-2 0.2797639049518 

-1 0.2897560966969 0.2897564200788 0.2897567439622 

0.2997483666777 0.2997486993479 0.2997490325375 0.2997493662436 0.2997497004638 

1 0.3097414024170 0.3097417452810 0.3097420886767 

2 0.3197345612229 



TABLE III: L on a regular (r* = 0.4 + ih, 0, = 0.3 + jh) grid with h = 0.01, m = 1 and a = 0.1; used in the metrics 11391 and 
PI. 

i -2 -1 1 2 

J 

-2 10.5859871980039 

-1 10.2472057969004 10.2608136646930 10.2744788974703 

9.9325339379316 9.945685737941 9.9588930170935 9.9721559726835 9.9854748024450 

1 9.6651644706586 9.6779990517326 9.6908877382673 

2 9.4161849439047 
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TABLE IV: 9H( ggf" ) = wb/3 on a regular (r* = 0.4 + ih, 0* = 0.3 + jh) grid with h = 0.01, m = 1 and a = 0.1; used in the 
metric 13911 . 

i -2 -1 1 2 

-2 0.0000209409007 

-1 0.0000216623673 0.0000216050272 0.0000215476741 

0.0000223786889 0.0000223196230 0.0000222605430 0.0000222014497 0.0000221423438 

1 0.0000229679828 0.0000229071865 0.0000228463764 

2 0.0000235446994 



TABLE V: The components of R a p for the metric (I39t at (r» = 0.4, 9, = 0.3) when m = 1, a = 0.1 for the 3 values of h shown. 
The data is fitted to R a p[h] — R a p[0] + h n B a p and the values of the 3 unknowns n (the convergence rate, shown below), 
Rap^t] (the- convergence limit, shown below), and B a p (not shown) can then be found. In the last column we assume that the 
convergence limit is zero, and the data from the two smallest values of h is fitted to R a p\H\ = h p C a p; then the convergence 
rate p is found and shown below. 



Raf3 

Roo 
Rm 

R 2 
R 03 

Rio 
Rn 

Rl2 
Rl3 

R20 
R21 
R22 
R23 
R30 
R31 
R32 
R33 



h = 0.01 
-.0000000380807 
-.0000000356994 
.0000000002564 
.0000001685141 
-.0000000356994 
.0000001077265 
-.0000000012202 
.0000000005647 
.0000000002564 
.0000003176420 
.0129878554048 
-.0000000137600 
.0000001685141 
.0000000005647 
-.0000000137600 
-.0000111845769 



h = 0.02 
.0000001511411 
.0000001416176 
.0000000010245 
.0000006739409 
.0000001416176 
.0000004299945 
.0000000044084 
.0000000022268 
.0000000010245 
.0000012706415 
.0521063939839 
.0000000550017 
.0000006739409 
.0000000022268 
.0000000550017 
.0000447327548 



h = 0.04 
.0000006048914 
.0000005668227 
.0000000040897 
.0000026940829 
.0000005668227 
.0000017202702 
.0000000175449 
.0000000089078 
.0000000040897 
.0000050765232 
.2109476828529 
.0000002195263 
.0000026940829 
.0000000089078 
.0000002195263 
.0001788197801 



n 

2.0048 
2.0052 
1.9966 
1.9989 
2.0052 
2.0013 
2.0428 
2.0071 
1.9966 
1.9977 
2.0217 
1.9961 
1.9989 
2.0071 
1.9961 
1.9989 



R a gM 
-.0000000005609 
-.0000000005629 
-.0000000000005 
-.0000000001356 
-.0000000005629 
.0000000004374 
-.0000000001985 
.0000000000143 
-.0000000000005 
-.0000000007055 
.0002061572158 
.0000000000366 
-.0000000001356 
.0000000000143 
.0000000000366 
.0000000099046 



P 



1. 

1.9880 
1.9985 
1.9998 
1.9880 
1.9969 
1.8531 
1.9794 
1.9985 
2.0001 
2.0043 
1.9990 
1.9998 
1.9794 
1.9990 
1.9998 



TABLE VI: The L 2 -norm \\R a p\\ of the metric at the gridpoints (r* = 0.4,0.5, 9, = 0.3,0.6,0.9,1.2) when m = 1, 

a — 0.1, 0.2 for the 3 values of h shown. The data is fitted to ||J? a| g[/i] || — \\R a p [0] || + h n B a p and the values of the 3 unknowns n 
(the convergence rate, shown below), ||_R a( 3[0]|| (the convergence limit, shown below), and B a p (not shown) can then be found. 
In the last column we assume that the convergence limit is zero, and the data from the two smallest values of h is fitted to 
R a p[h] — h p C a p; then the convergence rate p is found and shown below 



0.2 
0.2 
0.2 
0.2 
0.2 
0.2 
0.2 
0.2 



r* a, 

0.4 0.3 

0.4 0.6 

0.4 0.9 

0.4 1.2 

0.5 0.3 

0.5 0.6 

0.5 0.9 

0.5 1.2 

0.4 0.3 

0.4 0.6 

0.4 0.9 

0.4 1.2 

0.5 0.3 

0.5 0.6 

0.5 0.9 

0.5 1.2 



h = 0.01 
.0032469650568 
.0002130833763 
.0000389028104 
.0000293455458 
.0032470081286 
.0002131086500 
.0000389167189 
.0000293127899 
.0032453323914 
.0002143149911 
.0000393784101 
.0000298331343 
.0032454107976 
.0002143150084 
.0000393833905 
.0000297893272 



h = 0.02 
.0130266033030 
.0008527906803 
.0001556196190 
.0001173591202 
.0130267772029 
.0008528867283 
.0001556755506 
.0001172283812 
.0130201489364 
.0008577236053 
.0001575247969 
.0001193103237 
.0130204532491 
.0008577193834 
.0001575427252 
.0001191387256 



h = 0.04 
.0527369396876 
.0034184910131 
.0006226264251 
.0004691359914 
.0527376265713 
.0034188769168 
.0006228548385 
.0004686121003 
.0527123284154 
.0034382789181 
.0006302707322 
.0004769318246 
.0527135097953 
.0034382574117 
.0006303423768 
.0004762466973 



n 

2.0217 
2.0038 
2.0004 
1.9988 
2.0217 
2.0039 
2.0004 
1.9989 
2.0217 
2.0038 
2.0005 
1.9988 
2.0217 
2.0039 
2.0004 
1. 



[o] II 

.0000515392842 
.0000006083421 
.0000000127093 
-.0000000231698 
.0000515376608 
.0000006107702 
.0000000130975 
-.0000000234006 
.0000516338458 
.0000006116897 
.0000000140940 
-.0000000245470 
.0000516341443 
.0000006130839 
.0000000149722 
-.0000000260206 



P 

2.0043 
2.0008 
2.0001 
1.9997 
2.0043 
2.0008 
2.0001 
1.9997 
2.0043 
2.0008 
2.0001 
1.9997 
2.0043 
2.0008 
2.0001 
1.9998 
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TABLE VII: The L 2 -norm \\R a p\\ for the PI metric at the grid points (r, = 0.4,0.5, 0* = 0.3,0.6,0.9,1.2) when m = 1, 
a = 0.1, 0.2 for the 3 values of h shown. The data is fitted to ||7? Q| a[/i] | = \\R a g [0] 1 1 — I — h"B a/ 3 and the values of the 3 unknowns n 
(the convergence rate, shown below), \\R a p [0] || (the convergence limit, shown below), and B al 3 (not shown) can then be found. 
In the last column we assume that the convergence limit is zero, and the data from the two smallest values of h is fitted to 
R a p[h] — h p C a f}\ then the convergence rate p is found and shown below 



a r* v * 

0.1 0.4 0.3 

0.1 0.4 0.6 

0.1 0.4 0.9 

0.1 0.4 1.2 

0.1 0.5 0.3 

0.5 0.6 

0.5 0.9 

0.5 1.2 

0.2 0.4 0.3 

0.2 0.4 0.6 

0.2 0.4 0.9 

0.2 0.4 1.2 

0.2 0.5 0.3 

0.2 0.5 0.6 

0.2 0.5 0.9 

0.2 0.5 1.2 



0.1 
0.1 
0.1 



h = 0.01 
.0032469652386 
.0002130840309 
.0000389034943 
.0000293439533 
.0032470083035 
.0002131092800 
.0000389173773 
.0000293112574 
.0032453330980 
.0002143175794 
.0000393812217 
.0000298268109 
.0032454114777 
.0002143174996 
.0000393860967 
.0000297832468 



h = 0.02 
.0130266040297 
.0008527932981 
.0001556223556 
.0001173527506 
.0130267779027 
.0008528892486 
.0001556781866 
.0001172222518 
,0130201517634 
.0008577339580 
.0001575360447 
.0001192850311 
,0130204559705 
.0008577293475 
,0001575535526 
,0001191144036 



h = 0.04 
.0527369425946 
,0034185014848 
.0006226373828 
.0004691105137 
.0527376293706 
,0034188869986 
.0006228653942 
,0004685875840 
.0527123397244 
.0034383203299 
,0006303157634 
.0004768306609 
.0527135206818 
.0034382972686 
.0006303857252 
.0004761494153 



n 

2.0217 
2.0039 
2.0004 
1.9988 
2.0217 
2.0039 
2.0004 
1.9989 
2.0217 
2.0039 
2.0005 
1.9988 
2.0217 
2.0039 
2.0005 
1.9988 



\\R«p[o]\\ 


P 


.0000515392845 


2.0043 


.0000006083423 


2.0008 


.0000000127102 


2.0001 


-.0000000231700 


1.9997 


.0000515376608 


2.0043 


.0000006107701 


2.0008 


.0000000130977 


2.0001 


-.0000000234007 


1.9997 


.0000516338460 


2.0043 


.0000006116900 


2.0008 


.0000000140979 


2.0001 


-.0000000245469 


1.9997 


.0000516341442 


2.0043 


.0000006130844 


2.0008 


.0000000149752 


2.0001 


-.0000000260198 


1.9998 



TABLE VIII: The L 2 -norm \\R af} \\ for the FL metric at the grid points (r = 4, 5, = 0.3, 0.6, 0.9, 1.2) when m = 1, a = 0.1, 0.2 
for the 3 values of h shown. The data is fitted to ||.R a/ 3[/i]|| = ||i? Q ^[0]|| + h n B a p and the values of the 3 unknowns n (the 
convergence rate, shown below), ||-R a /3[0]|| (the convergence limit, shown below), and B a p (not shown) can then be found. 
In the last column we assume that the convergence limit is zero, and the data from the two smallest values of h is fitted to 
R a p[h] = h p C a f3; then the convergence rate p is found and shown below 



a 


r 


6 


h = 0.001 


h = 0.002 


h = 0.004 


n 


ll#a,3[0]|| 


P 


0.1 


4 


0.3 


.0000457585190 


.0001830406822 


.0007322693156 


2.0003 


.0000000089023 


2.0001 


0.1 


4 


0.6 


.0000025606471 


.0000102426033 


.0000409714183 


2.0000 


.0000000001052 


2.0000 


0.1 


4 


0.9 


.0000004382731 


.0000017531700 


.0000070127199 


2.0000 


-.0000000000300 


2.0001 


0.1 


4 


1.2 


.0000002993400 


.0000011973612 


.0000047894342 


2.0000 


-.0000000000017 


2.0000 


0.1 


5 


0.3 


.0000426162736 


.0001704709495 


.0006819789861 


2.0003 


.0000000079718 


2.0000 


0.1 


5 


0.6 


.0000024663873 


.0000098656545 


.0000394635674 


2.0000 


.0000000000586 


2.0000 


0.1 


5 


0.9 


.0000004274473 


.0000017097918 


.0000068391881 


2.0000 


.0000000000012 


2.0000 


0.1 


5 


1.2 


.0000002969713 


.0000011878357 


.0000047512974 


2.0000 


.0000000000169 


1.9999 


0.2 


4 


0.3 


.0000666096949 


.0002664508787 


.0010659963373 


2.0003 


.0000000160412 


2.0001 


0.2 


4 


0.6 


.0000030950589 


.0000123802850 


.0000495225364 


2.0000 


.0000000001332 


2.0000 


0.2 


4 


0.9 


.0000004998811 


.0000019995370 


.0000079982078 


2.0000 


.0000000000011 


2.0000 


0.2 


4 


1.2 


.0000003097359 


.0000012388778 


.0000049554632 


2.0000 


.0000000000239 


1.9999 


0.2 


5 


0.3 


.0000570798645 


.0002283289638 


.0009134670723 


2.0003 


.0000000125727 


2.0001 


0.2 


5 


0.6 


.0000028663076 


.0000114652931 


.0000458624019 


2.0000 


.0000000001087 


2.0000 


0.2 


5 


0.9 


.0000004733989 


.0000018936611 


.0000075746822 


2.0000 


-.0000000000248 


2.0000 


0.2 


5 


1.2 


.0000003046906 


.0000012187660 


.0000048750369 


2.0000 


-.0000000000045 


2.0000 



